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ABSTRACT 


The  purpose  of  this  research  is  to  improve  our  understanding 
of  shock  initiation  and  propagation  of  detonation  in  granular  or 
heterogeneous  solids.  Reactions  in  initiating  shocks  were  studied 
by  measuring  pressure  and  specific  volume  in  the  shocked  state  (2  to 
100  kbar)  for  the  physically  similar  granular  reactive  and  nonreactive 
aggregates  ammonium  perchlorate  (a  low-power,  insensitive  explosive) 
and  potassium  chloride  (an  inert  comparison  material) .  Longitudinal 
sound  velocities  in  compacts  of  the  same  material ,  over  a  range  of 
bulk  densities  and  grain  sizes,  were  also  measured.  Loci  of  shocked 
states  for  NI^CIO^  and  KC1  aggregates  show  a  dramatic  qualitative 
difference,  leading  us  to  the  conclusion  that  reaction  occurs  at  the 
shock  front.  This  conclusion  opens  the  possibilities  that  nonreac¬ 
tive  shocks  in  granular  explosives  and  propellants  do  not  exist,  that 
the  amount  of  initial  reaction  is  related  to  the  shock  pressure,  and 
that  mechanisms  of  buildup  to  steady  detonation  must  take  into 
account  reaction  from  time  of  shock  entry. 

Finally,  a  dependence  of  transition  behavior  on  grain  size  was 
demonstrated;  the  faster  transition  from  shock  to  detonation  occurs 
for  smaller  grain  size,  i.e.,  larger  grain  surface  per  unit  volume. 
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INTRODUCTION 


This  document  reports  work  carried  out  between  1  January  1966  and 
31  December  1966  on  Contract  AF  49(638) -1669 ,  "Detonation  Characteristics 
of  Low  Density  Granular  Materials." 

The  purpose  of  the  research  is  to  develop  information  leading  to  an 
understanding  of  initiation  and  propagation  of  detonation  in  granular 
materials  under  conditions  of  low  bulk  density.  Under  these  conditions 
the  specific  effects  related  to  heterogeneity  become  dominant  and  amenable 
to  investigation.  Effects  specific  to  heterogeneity  play  a  role  in  the 
behavior  of  most  practical  explosives  and  propellants.  Thus  the  informa¬ 
tion  obtained  will  be  useful  in  developing  a  knowledge  of  how  best  to 
handle  and  employ  solid  propellants  and  explosives.  Specifically,  it 
will:  give  insight  into  the  shock  sensitivity,  stability,  and  detonabil- 

ity  of  solid  p  'ipellants  and  explosives;  provide  background  information 
making  possible  the  design  of  more  reliable  igniters  and  detonators; 
aid  in  establishing  safe  manufacturing  procedures;  and  help  establish 
conditions  for  appropriate  propellant  and  explosive  use  in  aircraft, 
missiles,  and  space  ships. 

Work  this  year  has  been  concentrated  on  studying  the  mechanisms  by 
which  a  shock  may  initiate  a  detonation  in  a  granular  material  and  grow 
to  a  steady  detonation.  Toward  this  end,  two  related  experimental  pro¬ 
grams  have  been  undertaken.  In  the  first,  reactions  in  initiating  shocks 
have  been  studied  by  measuring  pressure  and  specific  volume  in  the  shocked 
state  for  the  physically  similar  granular  reactive  and  nonreactive  aggre¬ 
gates  ammonium  perchlorate  (a  low-power,  insensitive  explosive)  and 
potassium  chloride  (an  inert  comparison  material).  In  the  process  of 
carrying  out  these  experiments,  additional  data  on  the  shock-to-detonation 
transition  were  obtained  in  further  pursuance  of  an  effort  begun  in  the 
previous  year. 

In  the  second  program,  measurements  of  longitudinal  sound  velocity 
were  made  in  columns  of  granular  ammonium  perchlorate  and  potassium 
chloride  over  similar  ranges  of  density  and  grain  size.  The  data  are 
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needed  to  complete  the  shock  compression  curves  because  the  longitudinal 
sound  velocity  is  the  limiting  wave  velocity,  for  a  wave  of  infinitesimal 
pressure  and  material  velocity.  Measurements  were  made  over  a  bulk 
density  range  from  1.0  to  1.9  g/cc  and  for  three  grain  distributions, 
about  mean  diameters  of  13,  76,  and  500ji.  Also  measured  for  each 
material  was  the  longitudinal  sound  velocity  along  one  axis  of  a  single 
crystal . 

In  Part  II  we  report  the  work  on  the  longitudinal  sound  velc  y. 

In  Part  III  we  describe  the  shock  compression  work,  show  how  the  sound 
velocity  results  are  utilized  in  the  study  of  the  behavior  of  entering 
shocks,  and  interpret  the  results  in  terms  of  chemical  reaction  within 
the  shock  and  of  detonation  initiation.  In  addition,  the  dependence  of 
detonation  transition  behavior  on  grain  size  is  shown. 


II  LONGITUDINAL  SOUND  VELOCITIES  OF  GRANULAR  COMPACTS  OF  AMMONIUM 


PERCHLORATE  AND  POTASSIUM  CHLORIDE 

Measurements  were  made  of  longitudinal  sound  velocity  in  columns  of 
granular  ammonium  perchlorate  and  potassium  chloride  over  similar  ranges 
of  density  and  grain  size.  The  work  is  part  of  the  effort  to  determine 
the  locus  of  shocked  states  of  physically  similar  granular  reactive  and 
nonreactive  aggregates,  and  the  relationship  of  these  loci  to  initiation 
of  detonation  and  chemical  reaction  in  shocks  in  solid  propellants  and 
explosives.1  The  longitudinal  sound  velocity  is  the  limiting  wave  velocity, 
for  a  wave  of  infinitesimal  pressure  and  material  velocity.  The  ammonium 
perchlorate,  a  low-power  insensitive  explosive,  was  the  reactive  material; 
the  potassium  chloride  was  the  physically  similar  nonreactive  material 
studied  for  comparison.  Measurements  were  made  over  a  bulk  density  range 
from  1.0  to  1.9  g/cc  and  for  three  grain  distributions,  about  mean  diameters 
of  13,  76,  and  500u.  Also  reported  for  each  material  is  the  longitudinal 
sound  velocity  along  one  axis  of  a  single  crystal. 

The  velocity  is  calculated  from  a  determination  of  the  time  required 
for  an  acoustic  disturbance  to  traverse  a  known  thickness  of  sample,  a 
method  previously  described  by  Goettelman  and  Evans.2  Pulses  are  pro- 
duced  by  a  PZT4  transducer,  12  mm  in  diameter,  energized  by  400  volt 
square  electrical  pulses  of  0,6  usee  duration.  To  eliminate  complications 
associated  with  transducer-specimen  coupling,  several  thicknesses  of 
specimens  are  used  and  the  slope  of  a  least  squares  fit  straight  line 
through  the  sample  thickness-transit  time  points  is  taken  to  be  an 
estimate  of  the  infinite  wave  length  longitudinal  wave  velocity. 

The  materials  were  propellant  grade  ammonium  perchlorate  and  reagent 
grade  potassium  chloride.  The  material  was  ground  or  tumbled  to  produce 
approximately  spherical  particles,  after  which  it  was  dried  in  an  oven 
at  100°C  for  12  hours.  Particle  size  c  stribution  was  determined  with  a 


*  Lead-zirconate-titanate  piezoelectric  ceramic,  manufactured  by  Clevite 
Corporation,  Bedford,  Ohio. 


Sharpies  Micromerograph  (a  Stokes'  law  particle  size  analyzer).  Mean 
diameter  13/2  particles  had  a  size  distribution  ranging  from  3  to  100(2, 
with  50%  below  13 a.  Mean  diameter  76(2  particles  were  in  a  distribution 
ranging  from  53  to  104(2  with  50%  below  76(2.  Mean  diameter  500a  particles 
were  in  a  distribution  with  90%  at  500/2  and  the  remainder  in  coarser  and 
finer  sizes. 

Low  density  granular  compacts  at  1.00  and  1.20  g/cc  were  prepared 
by  packing  directly  in  a  styrofoam- 'Lucite '  tube  assembly  which  aligned 
the  transmitting  and  receiving  transducers.  Weighed  amounts  of  material 
were  placed  in  the  15.7  mm  I.D.  tube  and  the  upper  transducer  was  pressed 
down  to  accurately  measured  scratch  lines  to  arrive  at  the  desired  bulk 
density.  Errors  in  volume  and  weight  measurements  indicate  a  probable 
density  excursion  of  ±  0.05  g/cc.  Increments  of  length  were  obtained  by 
adding  increments  of  material,  and  alternatively  by  emptying  the  tube  and 
reloading  the  total  amount.  Comparison  runs  of  the  two  loading  methods 
showed  no  significant  difference  in  the  measured  transit  times.  A  velocity 
was  determined  from  transit  time  measurements  on  at  least  10  sample 
thicknesses  ranging  from  3  to  24  mm. 

High  density  cylindrical  pellets  of  4  mm  length  and  12.7  mm  diameter 
were  prepared  by  pressing  to  1.90  ±  0.02  g/cc  in  a  die.  The  pellets  were 
hard  and  translucent,  suggesting  that  the  original  grain  size  distributions 
were  altered  by  crushing.  A  velocity  was  determined  from  transit  time 
measurements  on  at  least  five  sample  thicknesses,  ranging  from  4  to  20  mm, 
obtainea  by  stacking  an  appropriate  number  of  pellets. 

Crystal  sound  velocities  were  measured  on  samples  consisting  of 
stacked  crystals.  The  KC1  crystals  (obtained  from  Harshaw  Chemical 
Company)  were  optical  quality  right  cylinders,  15  mm  dia.  x  3  mm  thick; 
the  axes  of  the  cylinders  were  in  the  [100]  direction,  and  the  measure¬ 
ment  was  made  in  the  axial  direction.  The  NH4CIO4  crystals  (generously 
furnished  by  Professor  E.  E.  Petersen  of  the  University  of  California 
at  Berkeley)  were  approximately  2  mm  thick  by  12  mm  wide.  Ammonium 
perchlorate  crystallizes  in  the  orthorhombic  form.  The  longitudinal 


wave  velocity  was  measured  perpendicular  to  a  broad  face,  in  the  direction 
of  tne  7.453  %.  spacing. **  A  thin  coating  of  petroleum  jelly  was  applied 
to  the  crystal  surfaces  to  improve  the  transmission  between  stacked 
crystals  when  they  were  not  perfectly  plane. 

Table  I  summarizes  results  in  terms  of  longitudinal  sound  velocities 
±  2  (7 ,  where  (7  is  the  standard  deviation  from  the  velocity  given  by  a 
straight  line  least  squares  fit.3  Figure  1  is  a  plot  of  sound  velocity 
versus  bulk  density.  Some  previously  published  data  points  are  included 
for  comparison.4 


Table  I 

LONGITUDINAL  SOUND  VELOCITIES  OF 
GRANULAR  COMPACTS  OF  NH4CIO4  AND  KC1 


Density 

Longitudinal 

Sound  Velocity  (mm/usec) 

NH4C104 

Mean  Particle 

Size 

KCl  Mean  Particle  Size 

(g/cc) 

500 fj, 

76  /I 

13  U 

76  /I 

134 

1.00 

* 

0 . 32*0 . 05 

0.72*0.03 

0.42*0.03 

0.53*0.04 

1.20 

0.57*0.08 

0.49*0.04 

1.21*0.13 

1.36*0.14 

1.18*0.08 

1.55 

1.79±0.36 

— 

2.26*0.12 

— 

— 

1.90 

2.18*0.56 

2.97*0.12 

3.30*0.06 

3.53*0.20 

3 . 87*0 . 33 

Crystal 

3.90*0.22  in  direction  of  7. 

0 

453  A  spacing 

4.70*0.34  in 

[lOO]  direction 

*  Below  pour  density. 


Longitudinal  sound  velocity  is  proportional  to  the  square  root  of 
a  stiffness  factor  (bulk  modulus  +  4/3  shear  modulus)  and  inversely 
proportional  to  the  square  root  of  an  inertia  factor  (bulk  density) .5 
While  models  of  aggregates  have  been  developed  to  estimate  the  stiffness 
factor  and  thereby  the  sound  velocity,6  they  require  knowledge  of 


**  Verified  with  X-ray  diffraction  by  Paul  S.  DeCarli . 
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properties  such  as  the  crystal  elastic  moduli  to  obtain  numerical 
results.  In  the  absence  of  this  information  foi  NH4CIO4  ,  we  discuss 
the  results  only  qualitatively. 

The  stiffness  and  hence  the  sound  velocity  of  a  porous  pressing 
would  be  expected  to  decrease  with  decreasing  bulk  density,  tending 
toward  zero  at  pour  density.  The  data  confirm  this  expectation.  The 
stiffness  of  a  pressing  might  be  expected  to  depend  not  so  much  on  par¬ 
ticle  size  as  on  size  distribution,  texture,  and  shape.  Sound  velocity 
in  the  KC1  pressings  is  found  to  be  independent  of  particle  size  within 
the  standard  deviation,  while  that  for  NH4CIO4  pressings  exhibits  a 
particle  size  dependence  at  each  bulk  density.  Even  at  1.9  g/cc  or 
97.5%  of  crystal  density,  where  materials  of  various  original  particle 
sizes  might  be  expected  to  have  nearly  the  same  stiffness  factor  and  thus 
nearly  the  same  sound  velocity,  the  particle  size  effect  persists.  Sub¬ 
stantial  increase  and  variability  in  sound  velocity  were  observed  for 
pressings  made  from  material  not  completely  dry;  the  moisture  caused 
caking,  which  increased  the  rigidity  of  the  aggregate. 


Ill  REACTION  IN  INITIATING  SHOCKS  IN  GRANULAR  AMMONIUM  PERCHLORATE 


A .  Shock  Compression  of  Reactive  and  Nonreactive  Granular  Materials 

As  part  of  our  study  of  shock  initiation  of  detonation  in  granular 
propellants  and  explosives,  pressure  and  specific  volume  in  the  shocked 
state  were  measured  for  the  physically  similar  granular  reactive  and 
nonreactive  aggregates  ammonium  perchlorate  and  potassium  chloride.  Shock 
compression  curves  are  reported  for  aggregates  of  the  two  materials  at 
1.0  cc/g  initial  specific  volume,  with  a  mean  particle  size  diameter  of 
13u.  Initial  shock  velocity  in  the  aggregate  and  free  surface  velocity 
of  the  driver  material  were  obtained  from  explosive-attenuator  and  light 
gas  gun  experiments.  The  data  extend  from  2.5  to  97.5  kbar. 

A  new,  explosive-driven,  multilayer  driver  system  was  designed  to 
provide  shocks  of  about  10  kbar  pressure  and  5  Usee  duration  in  the  ini¬ 
tial  density  1.0  g/cc  materials.  Data  were  also  obtained  using  a  similar 
type  of  driver  system,  developed  in  1965,  which  provides  shocks  of  2.5 
kbar;  these  experiments  were  performed  to  verify  the  low  pressure  region 
of  the  NH4C104  shock  compression  response  obtained  jn  the  previous  year. 

A  light  gas  gun  was  used  as  an  alternative  instrument  for  shock  experi¬ 
ments,  providing  an  independent  check  on  results  obtained  with  the 
explosive  driver  technique.  The  gun,  shown  in  Fig.  2,  uses  helium  to 
drive  aluminum  projectiles  at  velocities  up  to  700  m/sec  against  the 
material  being  tested.  A  projectile  is  shown  approaching  the  specimen¬ 
holding  assembly  at  the  end  of  the  gun  barrel.  The  projectile  motion  is 
monitoreu  by  pins  and  an  electronic  counter.  Arrival  at  the  face  of  the 
specimen  is  recorded  by  tilt  pins  which  lie  in  the  same  plane  as  the  face. 
A  shock  then  traverses  the  specimen  and  is  sensed  at  the  rear  face  by  a 
piezoelectric  or  resistance  wire  sensing  gage.  Special  electronic  cir¬ 
cuitry  is  used  to  measure  very  accurately  the  traverse  time,  which  is  the 
primary  measurement  required.  The  sensing  gage  response  also  tells  the 
shape  and  amplitude  of  the  shock  wave. 


■•jssjr’— r. 
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VACUUM  MANIFOLD 


GUN  BARREL 

GA-JBI9-4 


FIG.  2  CROSS  SECTION  OF  LIGHT  GAS  GUN  AT  IMPACT  END 


Results  expressed  in  terms  of  shock  pressure  and  final  volume  are 
shown  in  Fig.  3.  The  heavy  solid  line  is  for  NH4C104  at  initial  specific 
volume  V0  =  1  g/cc.  The  dashed  line  is  an  estimate  for  KC1  ^*ystai,  of 
initial  specific  volume  V  ,  extrapolated  from  several  points;  the  one 
indicated  was  obtained  in  this  laboratory.7  Open  points  indicate  data 
obtained  by  explosive  driver  techniques  and  solid  points  are  used  for 
data  from  gas  gun  techniques;  within  the  spread  of  experimental  accuracy, 
they  give  reasonable  confirmations.  For  KC1 ,  the  points  for  porous  and 
crystal  initial  states  lie  approximately  on  the  same  curve;  however,  the 
porous  material  starts  from  V0  and  must  approach  the  foot  of  the  curve 
at  some  small  slope.  The  single  point  for  porous  KC1  at  V0  =  1.25  cc/g, 
which  lies  to  the  right  of  the  crystal  curve,  is  from  Russian  experiments.8 

Preliminary  experiments  have  also  been  made  for  coarse  grain  size 
NH4CIO4  in  the  size  range  53  to  104u  diameter.  These  data  are  not 
reported  here,  pending  confirmation  and  analysis. 

The  data  of  Fig.  3  are  shown  in  Fig.  4  in  terms  of  shock  velocity 
versus  particle  velocity.  The  points  at  zero  particle  velocity  are  the 
longitudinal  sound  velocities  of  the  granular  materials  (see  Section  II). 


9 


VOLUME - cc/g 

GB-38I9-2B 

FIG.  3  SHOCK  COMPRESSION  OF'  POROUS  AND  SOLID  AMMONIUM 
PERCHLORATE  AND  POTASSIUM  CHLORIDE 
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PARTICLE  VELOCITY  —  mm/fisec 

8I-9SI9-IA 

FIG.  4  SHOCK  VELOCITY  VERSUS  PARTICLE  VELOCITY  FOR  SHOCK 
COMPRESSION  OF  POROUS  AMMONIUM  PERCHLORATE 
AND  POTASSIUM  CHLORIDE 


The  inert  KC1  points  lie  on  a  straight  line  extrapolating  to  the  sound 
velocity  point,  in  the  customary  behavior  of  inert  materials  and  solid 
explosives  at  pressures  below  those  at  which  reaction  occurs.9’10  A 
dashed  line  with  two  inflections  is  drawn  through  the  NH4CIO4  points. 

NH4CIO4  compacts  of  initial  specific  volume  V0  =  1.0  cc/g  show, 
at  a  computed  shock  pressure  of  2.5  kbar,  a  shocked  volume  which  is 
somewhat  larger  than  solid  crystal  volume  at  the  same  pressure;  they  also 
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show  a  positive  slope  for  dp/dv  between  2.5  and  10  kbar,  and  a  nega¬ 
tive  slope  at  higher  pressures.  The  computed  pressure  thus  is  double¬ 
valued  over  a  substantial  shocked  volume  range.  The  points  for  KC1  , 
at  grain  size  and  initial  bulk  density  similar  to  those  of  the  NH4CIO4 
samples,  show  much  smaller  volume  than  the  NH4C104  ,  and  the  points  are 
close  to  the  KC1  crystal  Hugoniot;  the  states  obtained  lie  on  a  curve 
going  up  and  to  the  left,  as  is  normal  for  inert  materials.  The  loci 
for  shocked  points  for  NH4CIO4  and  KC1  aggregates  thus  show  a  dramatic 
qualitative  difference.  V/e  interpret  this  difference  as  follows: 

The  data  points  for  the  compression  curves  were  determined  using 
entering  shock  velocities,  so  our  conclusions  are  concerned  with  the 
behavior  of  the  front  of  the  shock.  V/e  conclude  that  the  shock  front  in 
granular  reactive  material  does  not  resemble  a  shock  in  inert  material, 
even  at  shock  pressures  as  low  as  a  few  kilobars.  In  particular,  the 
reverse  curve  for  NH4CIO4  shown  in  Fig.  3  indicates  that  a  certain 
amount  of  reaction  occurs  at  the  shock  front.  These  results  open  the 
possibilities  that  nonreactive  shocks  in  granular  explosives  and  pro¬ 
pellants  do  not  exist,  that  the  amount  of  initial  reaction  is  related 
to  the  shock  pressure,  and  that  in  considering  the  mechanism  of  buildup 
to  steady  detonation  account  must  be  taken  of  chemical  reaction  from  the 
time  of  shock  entry  into  the  sample. 

B .  Shock-to-Detonation  Transition 

The  effect  of  grain  size  on  transition  to  detonation  was  explored 
in  NH4CIO4  by  using  a  grain  size  substantially  coarser  (53  to  104u 
diameter  range)  than  the  fine  grain  (13|i  mean  diameter)  previously  used. 
The  same  driver  system  previously  applied  to  the  fine-grained  material 
was  used  to  drive  a  plane-wave  shock  into  the  charge,  monitored  by  foil 
switches  in  one  case  and  by  fine  probe  wires  in  another.  Transition  to 
steady  velocity  appeared  at  50  to  75  mm  charge  length,  in  the  time  range 
20  to  26  flsec.  This  is  in  contrast  to  results  for  13U  material,  where 
the  transition  appears  at  about  15  mm  and  4  to  6  [X sec.1  This  comparative 
behavior  indicates  a  dependence  of  transition  uohavior  on  grain  size, 
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with  transition  occurring  more  quickly  with  smaller  grains.  One  expla¬ 
nation  of  this  result  would  be  that  transition  behavior  is  related  to 
grain  surface  area.  The  faster  trsnsition  with  smaller  grain  size  is 
then  explained  on  the  basis  of  larger  grain  surface  per  unit  volume. 
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